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We study the effects of thermal fluctuations on thin elastic filaments with sponta- 

neous curvature and torsion. We derive analytical expressions for the orientational 

correlation functions and for the persistence length of helices, and find that this 

length varies non-monotonically with the strength of thermal fluctuations. In the 

^r) ' weak fluctuation regime, the persistence length of a spontaneously twisted helix has 

(N ' 

JO ■ three resonance peaks as function of the twist rate. In the limit of strong fluctuations, 

■ all memory of the helical shape is lost. 

Recent advances in the art of micro-manipulation of molecules led to many experimental 
studies of the elasticity of biomolecules such as DNAihi, chromatin!, proteins!, and rodlike 
protein assembliesi~0. The tacit assumption behind many of the theories is that the elas- 
ticity of these biomolecules is of entropic originEHll and, consequently, they are modeled as 
random walks§. An alternative approach to the modeling of such systems is based on the 
assumption that the origin of elasticity is energy rather than entropy - there exists a lowest 
energy equilibrium configuration with associated spontaneous curvaturei , deviations from 
which give rise to restoring forces. While such an approach is a straightforward extension 
of the usual theory of elasticity of thin rodsR the description of arbitrary spontaneous cur- 
vature and twist involves rather complicated differential geometry and most DNA-related 
studies considered only fluctuations around the straight rod configuration!^ (see, however, 



references [18] and [19]). Following recent studies on the elasticity and stability of thin rods 
with arbitrary spontaneous curvature and torsionE3, in this work we investigate the effect 
of thermal fluctuations on the statistical properties of such filaments. We derive the differ- 
ential equations for the orientational correlation functions of the vectors pointing along the 
principal axes of the filament, and use them to calculate the correlators and the effective 
persistence length of an untwisted helix. Analytical expressions for the persistence length 
of a spontaneously twisted helix are obtained and it is found that this length varies non- 
monotonically with the amplitude of fluctuations and exhibits resonant-like dependence 
on the rate of twist. We would like to emphasize that although the present work is mo- 
tivated by recent studies of biomolecules, its aim is to construct a theoretical framework 
for the description of fluctuating string-like objects that goes beyond current models of 
polymer physics, and we do not attempt here to model particular experiments involving 
single-molecule manipulation. 

A filament of small but finite and, in general, non-circular cross-section, is modeled as 
an inextensible but deformable physical line parametrized by a contour length s (0 < s < L 
where L is the length of the filament). To each point s one attaches a triad of unit vectors 
t(s) whose component t 3 is the tangent vector to the curve at s, and the vectors ti(s) and 
t2(s) are directed along the axes of symmetry of the cross-section. Note that t(s), together 
with the inextensibility condition dx/ds = t 3 , gives a complete description of the space 
curve x(s), as well as of the twist of the cross-section about this curve. The rotation of the 
triad t as one moves along the curve is determined by the generalized Frenet equations 



Here flij(s) = J2k £ ijk^k(s), £ijk is the antisymmetric tensor and {ujk} are generalized tor- 
sions. 

The theory of elasticity of thin rod is based on the notion that there exists a stress- 
free reference configuration defined by the set of spontaneous curvatures and torsions {wofc}- 
The set {wofc} together with Eqs. (H) completely determines the equilibrium shape of the 
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filament. Neglecting excluded- volume effects and other non-elastic interactions, the elastic 
energy associated with some actual configuration {uJk} of the filament is a quadratic form 
in the deviations Sujk(s) = cJfe(s) — ujQk(s) 
kT r L 

U e i{Su k } = — dsJ2a k 5ul (2) 
2 Jo k 

where the bare persistence lengths a\ = EiIx/kT, a 2 = E\l 2 jkT and a 3 = E 2 (Ii + I 2 )/kT 
(T is the temperature and k is the Boltzmann constant) are expressed in terms of the elastic 
moduli Ei and the moments of inertia U about the axes of symmetry of the cross-sect ion0. 
The only limitation on the validity of Eq. (|2|) is that deformations are small on micro- 
scopic length scales, of order of the thickness of the filament. The elastic energy U e \ {Sook} 
determines the statistical weight of the configuration {uJk}- Calculating the corresponding 
Gaussian path integrals we find that {SuJi(s)) = and 

{Sui(s)Su)j(s')) = a^ 1 8 i j8(s — s') (3) 

We conclude that fluctuations of generalized torsions at two different points along the fila- 
ment contour are uncorrelated, and that the amplitude of fluctuations is inversely propor- 
tional to the corresponding bare persistence length. 

The statistical properties of fluctuating filaments are determined by the orientational 
correlation functions, (ti(s)tj(s')). In order to derive a differential equation for this correla- 
tors, we calculate the variation of t,(s) under the substitution s — > s + As. Integrating Eq. 
(H) yields in matrix notation (for small As): 

{/■s+As "| i-s+As i-s+As 

1 — J ds 1 n(s 1 ) + -J dsi J ds 2 0,(si)Q,(s 2 ) + 

X - £ +AS d Sl J** ds 2 [n(si)n(a 2 ) - ft(s 2 )0(s!)] + ... j t(s) (4) 

where the last term appears because of noncommutativity of matrices f2(si) and fl(s 2 ) for 
different si and s 2 . We multiply the above expression by tj(s'), average the result and 
note that for s + As > s > s' the fluctuations Sui(si) and 5ujj(s 2 ) are uncorrelated with 
the fluctuations of t»(s) and tj(s'). This implies that averages of products of fi's and t's 



factorize into products of the averages of JTs and those of t's. Since the averages of the 
terms in the square brackets in Eq. (|J) depend only on |si — S2I and their difference vanishes, 
in the limit As — ► this yields 

^ (uisfaW) = -£M«) MsfaW) (5) 

k 

where 

Kk = JiSik + £ iki UJ oi with 7i = H 7T T~~ ( 6 ) 

The above equations describe the fluctuations of filaments of arbitrary shape and flexibility 
and in the following this general formalism is applied to helical filaments. 

Consider a helix without spontaneous twist, such that the generalized spontaneous cur- 
vatures and torsions {c^ofc} are independent of position s along the contour. In order to 
visualize the stress-free configuration of such a filament, it is convenient to introduce the 
conventional Frenet triad of unit vectors which consists of the tangent, normal and binormal 
to the space curve spanned by the centerline, supplemented by a constant twist angle a 
which describes the rotation of the cross-section about this line. The relation between the 
two triads is given by u>o\ = Kocosao, ^02 = /tosina and = T o, where k and r are 
the constant curvature and torsion of the space curve. The rate of rotation of the centerline 
about the long axis of the helix is cuo = (kq + Tq) 1 ^ 2 , the helical pitch is 2titq/uJq and the 
radius is 2^^/^. For constant {kq, tq, ao}, A is a constant matrix and (tj(s)tj(s')) is given 
by ij—th element of the matrix exp [—A (s — s')]. The eigenvalues of the matrix A can be 
obtained analytically by solving for the roots of a characteristic cubic polynomial but the 
resulting expressions are cumbersome and will be presented in a longer report. Here we 
discuss only two limiting cases: 

Weak fluctuations, J2i 1% "C ojq In this case 

Us)t t (0)) = HK) e~ s/l + (l - ulJoot) cosMe^ 21 -^ (7) 

with the decay lengths I and 1$ defined by l^ 1 = J^klk^ok/^o anc ^ = Sfc a k 1 ° J okl ul o- 
The physical meaning of this correlator becomes clear by switching off thermal fluctuations 
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{ik = 0)- The first term on the rhs of this equation expresses the fact that the projection 
of any vector tj of the triad on the symmetry axis of the helix is constant, with magnitude 
ujQi/ujQ. The projections on the plane normal to this axis with magnitudes (1 — cu^/c^q) 1 / 2 , 
rotate with angular rate ujq. In the presence of weak fluctuations, the axis of symmetry of 
the helix becomes a random walk and the loss of correlations of its projections along the axes 
of the triad is described by the factor exp(— s/l). In the second term of Eq. (|7|), exp(— s/2l) 
describes the loss of correlations of the orientation of the plane normal to the axis of the helix. 
The angular persistence length 1$ results from averaging over the random angle of rotation 
((f)) with respect to the axis of the helix, (cosmos + <fi(s) — 0(0)]) = cos (a>os) exp (—s/21^) , 
with([0( S )-0(O)] 2 ) = S // . 

Strong fluctuations, ^ <^o- In this limit 

foOOt^O)) = e-^~<% (8) 

i.e., the correlators depend only on the bare persistence lengths and memory about the 
orientation of the vector t, decays over contour distance 7" 1 . Strong fluctuations destroy 
all phase coherence and all correlations between different vectors of the triad and lead to 
complete "melting" of the helix. 

We now proceed to calculate the effective persistence length l p which controls both the 
thermal fluctuations of a filament and its elastic response to external forces. It is defined as 
the ratio of the rms end-to-end separation (R 2 ) and the contour length of the filament L, 
in the limit L — > oo. The end-to-end vector is defined as R = J L t3 (s)ds and thus 

l p = lim | f L ds f ds' (t 3 (s)t 3 (s')) (9) 

L^oo L JO JO 

A simple calculation yields a result valid for arbitrarily strong fluctuations: 

l P = 2—s Tl± ^R (10) 

(71 cos 2 a + 72 sin a J + r 2 7 3 + 717273 

For non-vanishing curvature and torsion, this expression diverges in the weak fluctuation 
limit 7j — ► and the shape of the filament is nearly unaffected by fluctuations. Non- 
monotonic behavior is observed for "plate-like" helices, with large radius to pitch ratio, 



k /t . For 7, — > 0, the effective persistence length approaches zero (recall that L — > oo in Eq. 
(0)). Thermal fluctuations expand the helix by releasing stored length and initially increase 
the persistence length. Eventually, in the limit of strong fluctuations, the persistence length 
vanishes again (as 73" 1 ) because of complete randomization of the filament. The sensitivity 
to the constant angle of twist «o increases with radius to pitch ratio. In the opposite limit 
of "rod-like" helices Kq — > 0, the effective persistence length approaches 2/73 and becomes a 
function of a± and 02 only. Indeed, since straight inextensible rods do not have stored length, 
their end-to-end distance and persistence length are determined by random bending and 
torsional fluctuations only, and are independent of twist. 

The preceding analysis can be extended to fluctuating filaments with twisted stress-free 
states characterized by constant curvature k , torsion r and rate of twist of the cross-section 
about the centerline, dao/ds. The relation between generalized torsions {^ofc(-s)} and Frenet 
parameters {k , t , «o(s)} is given by u>oi = K o cos «o, 0J02 = sin a and u os = r + dao/ds. 
The calculation of the persistence length involves the solution of Eq. (|5|) with periodic 
coefficients. Details of the analytical solution will be given elsewhere. For filaments with 
circular cross-section a x = 02, the persistence length is independent of twist. In Fig. 1 we 
present a plot of the persistence length given in units of the helical pitch I* = I p Uq/2tyt , 
on the dimensionless rate of twist w* = 2ujQ 1 da>o/ds, for a "plate-like" helix with large 
radius to pitch ratio Ko/ T o and ribbon-like cross-section, a x <C 02- Curve 1 corresponds 
to the case of weak fluctuations, 7^ <C u . Throughout most of the range, the persistence 
length is independent of the rate of twist but a sharp peak appears at dao/ds = (see 
insert), accompanied by two smaller peaks at dao/ds = ±uj /2. Note that while in the limit 
of vanishing pitch, a ribbon-like untwisted helix degenerates into a normal ring, the cross- 
section of a twisted helix with dao/ds = ±u;o/2 rotates by ±7r and the helix degenerates into 
a Mobius ring. As the amplitude of fluctuations increases, the central peak transforms into 
a narrow dip and the two Mobius peaks become broad minima (curve 2). Further increase 
of fluctuations leads to the disappearance of the Mobius dips and the central dip becomes 
broad and shallow (curve 3). Finally when 7^ 3> uq, all dependence of the persistence length 

6 



on the spontaneous twist disappears (curve 4). Note that, as expected from the discussion 
following Eq. (PUD, the persistence length depends non-monotonically on the amplitude of 
thermal fluctuations. 

In order to understand the origin of the Mobius resonances we note that while the 
effective persistence length is a property of the space curve given by the Frenet triad, the 
microscopic Brownian motion of the filament arises as the result of random forces that act 
on its cross-section and therefore are given in the frame associated with the principal axes of 
the filament. Since the two frames are related by a rotation of the cross-section by an angle 
ceo(s), the random force in the Frenet frame is modulated by linear combination of sinao(s) 
and cosa ( s )- This gives a deterministic contribution to the persistence length which, to 
lowest order in the force, is proportional to the mean square amplitude of the random force 
and therefore varies sinosoidally with ±2«o(s). The observed resonances occur whenever 
the natural rate of rotation of the helix uq coincides with the rate of variation of this 
deterministic contribution of the random force, ±2dao/ds. 

In summary, we presented a statistical mechanical description of thermally fluctuating 
elastic filaments of arbitrary shape and flexibility. We would like to emphasize that the only 
limitation on the magnitude of fluctuations is that they are small on microscopic length 
scales, and that our model describes arbitrary deviations of a long filament from its equilib- 
rium shape. The general formalism was applied to helical filaments with and without twist. 
Strong thermal fluctuations lead to melting of the helix, accompanied by a complete loss of 
helical correlations. In general, the persistence length is a non-monotonic function of the 
amplitude of thermal fluctuations. Although through most of its range twist has a minor 
effect on the persistence length, resonant peaks and dips are observed whenever the rate 
of twist approaches zero or equals in absolute magnitude to half the rate of rotation of the 
helix. 
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Figure captions 

Plot of the dimensionless persistence length /* as a function of the dimensionless rate of 
twist w* for a helical filament with spontaneous curvature k = 1, and torsion r = 0.01 
(in arbitrary units). The different curves correspond to different bare persistence lengths: 
(1) ai = 100, a 2 = a 3 = 5000, (2) a x = 1, a 2 = a 3 = 100, (3) a x = 0.1, a 2 = a 3 = 10, (4) 
ai = 0.01, a 2 = a 3 = 10. A magnified view of the region of small twist rates is shown in the 
insert. 
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